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ABSTRACT 

The Darwin model of electromagnetic interaction is presented as 
a self-consistent theory, and is shown to be; an excellent approxima- 
tion to the Maxwell theory for slow electromagnetic waves. Since 
the fast waves of the Maxwell theory are absent, it is convenient 
for use in the computer simulation of the electromagnetic dynamics 
of nonrelativistic plasma. 
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Simulation studies of plasma behavior have usually been based 
either on the Coulomb model of particle interaction or on the full 
Maxwell theory. For the consideration of electromagnetic effects, it 
is not necessary to use the Maxwell theory, which includes the propa- 
gation of fast waves (whose group velocities are of the order of c ) , 
and which thus has the disadvantage of requiring a small time step. 
Indeed, the many plasma phenomena in which fast waves are unimportant 
are well described by the Darwin model^ of electrodynamics. 

The Darwin model, in which there is no retardation, has been used 

in the past to study electromagnetic interactions in microscopic 
2 

systems, and for the statistical mechanics of many-body systems, 

3 4 

both neutral and plasma. While for nonequilibrium plasma the Darwin 
model has been little used in analytic theoretical work, it seems most 
suitable for simulation studies, as pointed out by Hasagawa and Okuda,^ 
who rediscovered it in simulating one dimensional plasma dynamics 

In this paper we investigate the analytic properties of the Darwin 
model. The model is usually treated as a second order (in v/c) approxi- 
mation to the relativistic Maxwell theory.^ We present it here as a- 

self consistent theory arising from a particle Lagrangian: 

1 2 

L = S g nuv^ - C + M, where the kinetic energy is nonrelativistic, 
while the Coulomb energy, C = i Z e.0 5- + Z e.0 ex (r,t), and the magnetic 
energy M = | Z e i (v i /c)*A i + Z e ± (^,/e) *A ex (r ,t ), include both pair 
interactions and interactions with external potentials. The internal 
potentials are defined as 

Ar.Urp) h £ V r U ’ 

<j(A) 


(la) 
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- Y, (ib) 

j( 0 i) 

The magnetic interaction tensor is T(r) = i (l + rr )/r. The internal 
potential (lh) is solenoidal "because (d/3r)*T(r) = 0. We therefore 
choose the transverse or Coulomb gauge also for the external potential: 
V*A eX = 0. 

The canonical momentum of a particle 


£i 


dL/cfrv. 


m.v. + 




( 2 ) 


is a function of the velocities of all the particles through A" 1 '. The 
Lagrangian equation ~ Sh/chm, is the standard Newton-Lorentz equa- 
tion of motion, when the electromagnetic field is expressed in terms 
of the potentials. 

When the external potentials are static, the' invariance of L 

under time-translation implies conservation of energy dH/dt = 0, where 

1 2 

the Hamiltonian H = £ p. *v. - L has a simple form H = £ 77 m.v, + C + M 

’'■-i ~l d 11 

3 4. 

in terms of velocities (but not in terns of momenta 1 3 ). When the 
external field vanishes, the invariance of L under space -translation 
implies conservation of total canonical momentum: (d/dt)Z g.(t) = 0 . 

If the set of discrete particles is approximated by a continuum 

with charge and current densities p(r,t) and j(r,t), the internal 

i i in in 

potentials (0 ,A ) become fields 0 (r,t), A (r,t) satisfying 


- ^ 0 ln = kvp , 


(3a) 

- ^A ln = 4 tt^ t /c , 

V-A in = 0 . 

(3b) 



-4- 


UCRL-19859 


The som’ce of A ln is the transverse (or solenoidal) part of j, defined 
most simply in' terms of spatial Fourier transforms: 

^(k^t) = (j - kk)*j(k,t). Note that the retardationless elliptic 
equation (3b) has the same form as the Poisson equation (3a ) , whose 
solution is standard in simulation studies. 

Of the four Darwin equations for the internal electromagnetic 
field, only one: 

cV x B in = 4 tt^ + dE L ln /dt (4) 

differs from Maxwell's: the others are the same. Here E T = - is 

~L 

the longitudinal (or irrotat tonal) part of the electric field. The 

in 

longitudinal part of (4), 0 = 47Tj + Se /cit, expresses charge con- 

~L ~L 

servation in agreement with (3a). The transverse part of (4), 
cV x B ln = 4'ir^,, differs from Maxwell theory in omitting dE^'Vdt, 
and leads to the Biot-Savart law: B in (r, t) = c Vd 3 r'(v|r - r' | ^) x 

The differential conservation laws differ from those of the 
Maxwell theory. If we require the external field also to satisfy the 
Darwin field equations, the energy equation is 

-E-r = (S/St)[(|E L | 2 + B 2 )/8tt + Z s /d 3 v | m s v 2 f s ] + V-[(c/4tt)E x B 

+ (47rc)~ 1 (cA/at)(S0/at) + Z /d 5 v - m v 2 vf 1 , 

^ s 2 s ^ s J 

where f (r,v,t) is the particle distribution function, satisfying the 

S ^ 

standard Vlasov equation. The momentum equation is, in the absence 
of external fields, (3/8t)(c ^pA + Z g /d^v m vf g ) '= - V*P, where the 

■2 2 2 

momentum flux density is P s Z / d^v m wf h- l(B + |e t I + 2 A*fe r )/ 87 r 

~ s s~~ s =3 ~L ~ ~L 
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- f BB + E r E t + c -1 (Ai;. + E T A) Air. (Here E T h dE T /dt.) 

L — ~L-'L '~-L ~L~ J ~L ~L 

To see how the Darwin model modifies the plasma dynamics,, we con- 
sider linear oscillations of a uniform plasma. With ^(k,^) = 
n(k,ou) •E(k,u>) and the Darwin field equations, we obtain 


[e(k,co) - (n 2 + ljjj *E(k,a)) = 0 , 


(5) 


where e = I + "^g, rf s k^c^/cn^, J = I - kk, and a must be derived 

from the linearized kinetic equation. The dispersion equation o>(k'' 
results from the vanishing of the determinant of (?)• The Maxwell 
equations, on the other hand, yield 


2 2/2 


A A 


[e(k,o>) - n 2 ^]-E(k,<D) = 0. (6) 

For a cold plasma in a static magnetic field, e is independent of 
k. Then a simple relation 

n D 2 (co,e) = ii m 2 (cd,0) - 1 (7) 

exists between the Darwin and Maxwell refractive indices n^ and 

(0 is the angle between k and B .) Curves of n^ (cjq), for 0 = 0, 

0 = Tr/2, and 0 < 0 < 7r/2, may be found, e.g., in Refs. 7 and 8. From- 

these diagrams (and from the associated analysis), we see that for 

frequencies less than the highest resonant (n = oo ) frequency, \ n^ | 

js generally large compared to one, i.e., the phase velocity is much 

less than c. The exceptions are the small frequency intervals near 

the cutoffs (n~ = 0), where |n | < 1. We conclude that, except near 

the cutoffs, the Darwin model agrees well with the Maxwell results. 

2 

For frequencies above the highest resonance, n^ <1 and therefore 




/ 
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2 

rip <0. The supraluminous waves of the Maxwell theory, which are 

troublesome in plasma simulation, are replaced by evanescent waves in 

the Darwin model. In terms of an cd vs k diagram, the fast branch 

(oj -* kc as 1; -* oo ) simply disappears. 

For a hot plasma, where _e depends also on k, no simple relation 

like (7) can be found. However, comparison of (6) with (5) again leads 

2 

us to conclude that the Darwin model is reliable for slow waves (n » l) 

2 

and unreliable for fast waves (n < l). 

Finally, for an inhomogeneous plasma, the new modes (drift waves) 
have still lower frequencies and phase velocities, and hence the Darwin 
model should be valid here as well. 

We conclude that, in plasma simulation work, it is not necessary 
to choose between a Coulomb model and the full Maxwell theory; the 
Darwin model presents a third intermediate possibility. It describes 
electric and magnetic forces, includes induced fields as well as static 
fields, and yet retains much of the simple structure of the Coulomb 
model. Since the electrodynamics of a nonrelativistic plasma usually 
deals with slow waves, it is seen that the Darwin model has a wide 
range of applicability. 

In this discussion, we have made no attempt to deal with the 
formulation of differencing schemes, nor to face the difficulties 
pointed out by Hasegawa and Okuda concerning boundary conditions.^ 

We have benefitted greatly from the advice of C. K. Bird sail, 

J. Byers, A. Hasegawa, J. Killeen, and A. B. Langdon. 
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